Local density of states of ID Mott insulators and CDW states with a boundary 
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We determine the local density of states (LDOS) of one-dimensional incommensurate charge 
density wave (CDW) states in the presence of a strong impurity potential, which is modeled by a 
boundary. We find that the CDW gets pinned at the impurity, which results in a singularity in 
the Fourier transform of the LDOS at momentum 2kF. At energies above the spin gap we observe 
dispersing features associated with the spin and charge degrees of freedom respectively. In the 
presence of an impurity magnetic field we observe the formation of a bound state localized at the 
impurity. All of our results carry over to the case of ID Mott insulators by exchanging the roles of 
spin and charge degrees of freedom. We discuss the implications of our result for scanning tunneling 
microscopy experiments on spin-gap systems such as two-leg ladder cuprates. 
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In recent years scanning tunneling microscopy (STM) 
and spectroscopy (STS) techniques have proved to be 
very useful tools for studying strongly correlated elec- 
tron systems such as high temperature superconductors 
(HTSC)Jl|, carbon nanotubcs [2, and rarc-carth com- 
pounds 01 ■ The usage of magnetic tips has also enabled 
the investigation of magnetic properties . These exper- 
iments have motivated an intense research effort on the 
theory of STS the main focus being on quasiparti- 

cle properties in HTSC. STM measures the local differ- 
ential conductance dI/dV{x), which is directly related 
to the LDOS. Impurities break translational invariance 
and lead to a modification of the LDOS in their vicinity, 
from which one can infer characteristic properties of the 
bulk state of matter as well as the nature of its electronic 
excitations. Interestingly this holds even in strongly cor- 
related systems without quasiparticle excitations as was 
demonstrated for a non-magnetic impurity placed in a 
gapless Luttinger liquid [1, 0] , which can be regarded as a 
quantum critical CDW state. In this case, at low energies 
the impurity behaves effectively as a physical boundary 
Q where the phase of the CDW order parameter gets 
pinned, giving rise to induced CDW order. 

Here we consider STS in a ID strongly correlated sys- 
tem with a spin gap in the presence of an impurity. This 
problem is of interest to the study of (quasi) ID CDW 
systems, two-leg ladder materials with strong supercon- 
ducting correlations, and the stripe phases of HTSC @. 
We focus on the regime in which the scattering at the 
impurity is strong and hence at sufficiently low energies 
it acts as a physical boundary [l^l- The inherently non- 
perturbative nature of STS in gapped (quasi)-lD systems 
requires an entirely different treatment compared to pre- 
viously studied cases. Following Ref. Q we consider the 
spatial Fourier transform of the LDOS throughout, as 
this allows physical properties to be more easily identi- 
fied and is commonly used in experiments [l|, [1] . 

The starting point of our analysis is the continuum de- 



scription of a CDW state. The latter arises in two rather 
different kinds of ID correlated electron systems: (1) A 
partially filled band of spinful electrons coupled to optical 
phonons of frequency Wph- At energies small compared 
to the electron-phonon coupling results in an attrac- 
tive interaction between electrons that can overcome the 
Coulomb repulsion [ll| . (2)^rongly correlated two- and 
three- leg ladder systems [3, ■ Here, in spite of strongly 
repulsive electron interactions, a Mott state with a finite 
spin gap occurs for a range of dopings around half fill- 
ing [l3|. In both cases there is a broad range of param- 
eters such that at low energies the system gives rise to 
a CDW state characterized by a gapped spin sector and 
a gapless charge sector. Regardless of the microscopic 
origin of the spin gap, by taking the continuum limit and 
bosonizing one arrives at a spin-charge separated theory 
describing collective charge and spin degrees of freedom. 
We now imagine a potential impurity to be present which 
at low energies and temperatures effectively cuts the sys- 
tem into two disconnected parts [l^l • We then can model 
the impurity by a boundary condition on the continuum 
electron field ^'^(s; = 0) = 0, resulting in a spin-charge 
separated Hamiltonian of the form H = He + Hs 
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Here "I>c.s arc canonical Bose fields which satisfy the 
boundary conditions ^c,s{x = 0) = 0, and Oc,s are their 
dual fields. The charge and spin velocities fc,s, the Lut- 
tinger parameters Kc,s and coupling constant gs are func- 
tions of the parameters defining the underlying micro- 
scopic model. In the case of single chain electron-phonon 
systems we have Vg > Vc and > 1 as we are deal- 



ing with attractive electron-electron interactions. On the 
other hand, in the case of two-leg ladders there is a spin 
gap with Vs < Vc and A'c < 1. We will consider both sit- 
uations in what follows. The charge sector ([1]) describes 
gapless collective charge excitations propagating with ve- 
locity I'c, whereas the spin excitations are described by a 
sine-Gordon model with a boundary In the regime 
considered here the bulk spectrum of the latter consists 
of gapped (anti)soliton excitations. At the Luther-Emery 
point Ks = I/V2 the spin sector is equivalent to a free 
massive Dirac fermion [1^. As is well known, (IT|)-(l2|) re- 
duces to the low-energy theory of a half-filled one-band 
Mott insulator, provided we interchange charge and spin 
sector and then set Kg = 1 and kp = tt/2. By virtue of 
this connection all our results for CDW states carry over 
to ID Mott insulators. 

The central object of our study is the time-ordered 
Green's function in Euclidean space, 

(r, XUX2) = - (Obi Tr ^a{T, xi) (0, X2) |0b) . 

(3) 

Here |0b) is the ground state and t ~ it denotes imagi- 
nary time. At low energies the electron annihilation oper- 
ator can be decomposed into right- and left-moving com- 
ponents as ^aix) = e''''"^Ra{x) + e"""'^io-(a;), which 
reduces @ to 
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where e.g. G^,^, (Ob | i?,(T, xi) 4, (0, xa) | Ob). As 
we are interested in the LDOS, we ultimately want to set 
xi = X2- As was noted in Ref. @ it is useful to consider 
the Fourier transform of the LDOS as physical proper- 
ties can be more easily identified. In momentum space 
the RL and LR contributions occur in a different region 
{Q ~ ±2kF) compared to the RR and LL parts (Q ~ 0). 
In absence of a boundary we have G^^, ~ G^a' = as the 
charge parts of these Green's functions vanish. In pres- 
ence of a boundary left and right sectors are coupled and 
the Fourier transform of the Green's function ^ con- 
comitantly acquires a nonzero component at Q ~ ±2kp, 
which gives a particularly clean way of investigating im- 
purity effects. For this reason we focus on the 2fcF-part 
of the Green's function in what follows but note that the 
small momentum regime can be analyzed analogously. 

The Green's function G^^, factorizes into a product of 
correlation functions in the spin and charge sectors. The 
charge part can be determined by a standard mode ex- 
pansion [1^. The sine-Gordon model on the half- line ^ 
is known to be integrable for quite general boundary con- 
ditions This enables us to calculate the correlation 
functions in the spin sector using the boundary state for- 
malism introduced by Ghoshal and Zamolodchikov 
together with a form- factor expansion [l^. Taking into 
account only the leading terms wc arrive at (r > 0) 
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where gc^siT,x) are the contributions of the charge and 
spin sectors respectively. Here Kq is a modified Bessel 
function, the normalization constant Zi was obtained in 
Ref. [3 J and explicit expressions for the boundary reflec- 
tion amplitude K{9) are given in Ref. J/7| (at the Luther- 
Emery point wc have [I^l K{6) = i tanh|). The expo- 
nents in the charge sector arc related to the Luttinger 
parameters by a = (A'c + l/Kc)^/8, b = {K^ - l/A'c)^/8, 
and c = ~ Kf)/S. The subleading terms in §^ 

involve three or more particles in the intermediate state 
or higher orders in the boundary A'- matrix (20j . The 
Fourier transform of the LDOS for £' > is 

doc ^e'iEt^Q-)c{t,x,x), 

(7) 

where the Green's function has been analytically contin- 
ued to real time. For Q « 2k-p only contributes and 
using ([6]) we arrive at our main result 
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Here \q\ <^ 2kp, A = arcosh(^), Fi denotes Appell's 
hypergeometric function, hi{6) = 1, h2{9) = K(^9 + 
i f ) 6^/2, ui = 2(£; - Acosh6')/wcg -f isgn(wsg/A) S, and 
U2 = 2us(£'— Acosh6')/uc(ws'7— 2Asinh6')^-isgn(^;s(7/A — 
2 sinhfl) 6, where 6 0+. Below we plot for two dif- 
ferent parameter regimes. We smoothen the singularities 
in the LDOS by taking S small but finite. In exper- 
iments the singularities are broadened by instrumental 
resolution and temperature. The results presented below 
apply to the regime T <C E, A,Vc/ao (oq is the lattice 
spacing), where temperature effects are negligible. 

Repulsive Case: We first consider the case of repul- 
sive electron interactions {vg < Vc, < 1). In Fig. [1] 
we plot N„{E,2kp + q) for the case of unbroken spin 
rotational symmetry (ATg = 1). The LDOS is domi- 
nated by the strong peak at momentum 2kp, which has 
its origin in Ni{E,Q) and is indicative of the CDW or- 
der being pinned at the boundary. This is analogous to 
the Luttinger liquid case Q. At low energies above the 
spin gap A we further observe two dispersing features 
associated with the collective spin and charge degrees 



2 



E/A 



E/A 



FIG. 1: \N„{E,2kf + q)\ (arbitrary units) for = 0.8, 
Ks = 1 and Vc — 2vs. The curves are constant g-scans 
which have been offset along the y-axis by a constant with 
respect to one anothe r. We observe dispersing features at 
= + {vsq/2y and Ec = Vc\q\/2 + A. For g < the 

dispersing features are strongly suppressed. 



of freedom respectively: (1) a "charge peak" that fol- 
lows Ec = Vc\q\/2 + A and (2) a "spin peak" at position 
Eg = ■\/A2 + {vsq/2y. The charge peak arises from the 
contribution Ni{E,Q) to the Fomier transform of the 
LDOS, whereas the spin peak has its origin in N2{E, Q), 
which encodes the effects of the boundary on the spin de- 
grees of freedom. We note that like in the Luttinger liquid 
case the phase of exhibits characteristic jumps at 
the peak positions. 

Attractive Case: In Fig. [2] we show N„{E,2kY + q) 
for the case of attractive electron interactions {vs > Vc, 
Kc > 1) and unbroken spin rotational symmetry. The 
peak at 2kp is much less pronounced than in the repulsive 
case. We again observe charge and spin peaks that follow 
Ec and Eg respectively. For momenta q above a criti- 
cal value qQ = 2Avc/vsy^v^ — a third dispersing low- 
energy peak appears at E = Vc\q\/2 + A^/l — {vc/vg)^- 
This feature can be thought of as arising from a spin ex- 
citation with momentum go and a charge excitation with 
momentum g — go- This is reminiscent of what is found 
for the single-particle spectral function in the bulk [2]| . 

So far we have considered the simplest possible bound- 
ary conditions corresponding to a spin-independent phase 
shift TT. In reality one may expect an impurity to give 
rise to a local potential or local magnetic field, which 
result in more general phase shifts upon reflection of 
particles at the boundary. In particular, these more 
general boundary conditions can give rise to boundary 
bound states, see e.g. [l^. In order to exhibit the 
signature of boundary bound states in LDOS measure- 
ments we now consider boundary conditions of the form 
$c(r,0) = 0, <^s{t,0) = $° with < $° < tt. In terms 
of the right- and left-moving fields this corresponds to 
i?,(T,0) = -e-'/'*°/2L,(T,0), where /t = 1 = -fi- 
We note that these boundary conditions break spin ro- 



FIG. 2: \N^{E, 2k-e+q)\ (arbitrary units) for = 1.2, = 
1 and Vs — 2vc. W e observe dis persing features at Es, Ec, and 
E = Vc\q\/2 + Ay^l-ivc/v^y (for |g| > go only). 



tational symmetry. If we go over to the case of the ID 
Mott insulator by exchanging spin and charge degrees of 
freedom the spin rotational symmetry remains intact and 
the boundary conditions correspond to a local potential. 
A discussion of the general case with nontrivial boundary 
phase shifts in both spin and charge sectors is left for a 
longer publication. 

The leading terms in the form-factor expansion for 
the chiral Green's function G^^, are still given by ([6]), 
but now the boundary reflection amplitude depends on 
cr. For example, at the Luther-Emery point we have 
K^ie) = sin(i6i/2 - /^$o/2)/cos(i6i/2 + /a$°/2). We 
note that nevertheless the Green's function is still diag- 
onal in spin indices G^^, cx 5^^>. If 7r/2 < K^{9) has 
a pole in the physical strip < 3m 6' < tt/2, which gives 
rise to an additional term linear in the boundary reflec- 
tion matrix in the form-factor expansion. The resulting 
contribution to Ni{E,Q) has a non-dispersing singular- 
ity at the lower threshold A sin $g , see Fig. [31 The emer- 
gence of a non-dispersing feature within the gap signals 
the presence of a boundary bound state. We note that 
the boundary bound state appears only in the down-spin 
channel Ni{E, Q). On the other hand, if we were to con- 
sider N„{E < 0,Q), the additional feature would appear 
in the up-spin channel only. 

Finally we wish to discuss the implications of our re- 
sults for STM experiments on quasi-lD materials. Our 
results apply for energies above the 1D-3D cross-over 
scale, which is set by the strength of the 3D cou- 
plings. Perhaps the most interesting materials to which 
our findings may be applied are two-leg ladders like 
Sri4Cu2404i [1^. The model we have studied captures 
the essential features of the low-energy description of 
(weakly doped) two- leg ladders, namely a gapless charge 
sector and a gapped spin sector. The tunneling cur- 
rent measured in STM experiments is directly related 
to the local density of states N„{E,Q). From the 2kp- 
component of the current one can hence extract informa- 



tions about the CDW correlations in the presence of an 
impurity. In particular, we expect the various peaks in 
Na{E,Q) corresponding to the pinned CDW order and 
the dispersing spin and charge degrees of freedom to ap- 
pear in the Fourier transform of the local differential con- 
ductance. A possible asymmetry N-\{E,Q) — Ni{E,Q) 
may be detected using a magnetic STM tip 

In summary, we have determined the low energy LDOS 
in strongly correlated gapped ID systems such as Mott 
insulators and CDW states in the presence of a strong im- 
purity potential. We have shown that the spatial Fourier 
transform of the LDOS can be used to infer character- 
istic properties of the bulk state of matter. The LDOS 
is dominated by a singularity at 2k-p, which is indicative 
of the pinning of the CDW order at the position of the 
impurity. The LDOS further exhibits clear signatures 
of propagating collective spin and charge modes, which 
reflect the nature of the underlying electron-electron in- 
teractions. We have investigated the modification of the 
LDOS in the presence of impurity bound states and dis- 
cussed the relevance of our results to STM measurements 
on two-leg ladder materials like Sri4Cu2404i. 
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FIG. 3: |iV|(E,2fcF + q)\ (full line) and \Ni{E,2kv + q)\ 
(dashed line) for Vaq/A = 6, Kc = 1, Kb = l/\/2, Vb — 2vc, 
and $B = 0.9 TT. The singularity of Ni{E,2kF + q) at 
E = A sin $g is due to the boundary bound state. The broad 
maximum of at E ~ 1.8 is due to the excitation of the 
boundary bound state and additional charge excitations. 
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